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Interpretation and implication of the non-detection of GeV 
spectrum excess by Fermi 7-ray Space Telescope in most GRBs 
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ABSTRACT 

Since the launch of the Fermi Gamma-ray Space Telescope on 1 1 June 2008, significant de- 
tections of high energy emission have been reported only in six Gamma-ray Bursts (GRBs) 
until now. In this work we show that the lack of detection of a GeV spectrum excess in almost 
all GRBs, though somewhat surprisingly, can be well understood within the standard internal 
shock model and several alternatives like the photosphere-internal shock (gradual magnetic 
dissipation) model and the magnetized internal shock model. The delay of the arrival of the 
> 100 MeV photons from some Fermi bursts can be interpreted too. We then show that with 
the polarimetry of prompt emission these models may be distinguishable. In the magnetized 
internal shock model, high linear polarization level should be typical. In the standard inter- 
nal shock model, high linear polarization level is still possible but much less frequent. In 
the photosphere-internal shock model, the linear polarization degree is expected to be roughly 
anti-correlated with the weight of the photosphere/thermal component, which may be a unique 
signature of this kind of model. We also briefly discuss the implication of the current Fermi 
GRB data on the detection prospect of the prompt PeV neutrinos. The influences of the intrin- 
sic proton spectrum and the enhancement of the neutrino number at some specific energies, 
due to the cooling of pions (muons), are outlined. 

Key words: gamma rays: bursts — polarization — radiation mechanism: nonthermal — ac- 
celeration of particles — elementary particles: neutrino 



1 INTRODUCTION 

Gamma-ray Bursts (GRBs) are the most extreme explosion dis- 
covered so far in the universe. With the discovery of the after- 
glows and then the m easurement of the redshifts in 1997 (see 
Ivan Paradiis et al.l l200Ct for a review), the cosmological origin 
of GRBs has been firmly established. The modeling of the late 
(t > 10 4 s) afterglow data favors the external forward shock 
model (see |Piranll999l : lMeszaros l2002l ; IZhang & Meszaros |2004 
for reviews). The radiation mechanisms employed in the model- 
ing are synchrotron radiation and synchrotron self-Compton (SSC) 
scattering. In the early time the prolonged activity of the central 
engine plays an important role in pr oducing afterglow emission 
too, particularly in X-ray band (e.gjKatz et al. |[l998l : lFan & Weil 
l2005fclNousek et al. l2006l : IZhang et al. l2006h . The radiation mech- 
anisms, remaining unclear, are assumed to be the same as those 
of the prompt soft gamma-ray emission. It is expected that in the 
Fermi era the origin of the prompt emission can be better un- 
derstood. This is because the Large Area Telescope (LAT) and 
the Gamma-ray Burst Monitor (GBM) onboard Fermi satellite 
(http://fermi.gsfc.nasa.gOv/l can measure the spectrum in a very 
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wide energy band (from 8 keV to more than 300 GeV), with which 
some models may be well distinguished. For example, in the stan- 
dard internal shock model the SSC radiation can give rise to a dis- 
tinct GeV excess while in the magnetized outflow model no GeV 
excess is expected. 

Motivated by the detection of some > 100 MeV pho- 
tons from quite a few GRBs by the Com pton Gamma Ray 



Observatory satellite in 1991-2000 (e.g., iHurlev et 



mma Kay 
1I| [1994; 



iFishman & Meegan II 19951 : iGonzalez et al.'feo03h . the prompt high 
energy emission has been extensively investigated and most cal- 
culat ions are within the framework of the standard internal shocks 
(e.g. IPilla & Loeb ll998l:|Pe'er & Waxman l2004l : lGupta & Zhang] 



l2007l : lBosniaketal.ll2009l, cf. Giannios 2007). The detec tion 
prospect for LAT seems very promising (see lFan & Piranl2008l , for 
a recent review). 

Since the launch of Fermi satellite on 11 June 2008, sig- 
nificant detections of prompt high energy emission from GRBs 
have been only reported in GRB 080825C dBouvier et al]|2008[) . 
GRB 0809 16C |Abdo et al.l2009h. G RB 081024B JOmodei l2008h, 
GRB 090323 Johno et al. I l2009bh. GRB 090328 dCutini et al. I 
120091) . possibly and GRB 090217 dOhno et al. Il2009ah until now 
(5 May 2009). Though the detection of 3 prompt photons above 10 
GeV from GRB 080916C at redshift z ~ 4.5 dAbdo et alj|2009l : 
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iGreiner et al. 112009b is amazing and may imply a very high ini- 
tial Lorentz factor of the outflow Yi > 180 and an efficien t ac- 
celeration of particles to very high energy dZou et aLll2009r) . the 
non-detection of a significant > 100 MeV emission from most 
GRBs may be a better clue of the underlying physics. A delay 
in the onset of the > 100 MeV emission with respect to the soft 
gamma-rays, as detected in GRB 080825C, GRB 080916C, GRB 
081024B and GRB 090323, may be the oth er clue of the GRB 
physics dAbdo et all200^ : IOhno et al. 2009b). In this work we fo- 
cus on these two novel observational features. 

This work is structured as the following. In section 2 we dis- 
cuss the constraint of the current Fermi GRB data on the standard 
internal shock model and several alternatives. In section 3 we look 
for distinguished signals in linear polarization of the prompt emis- 
sion. In section 4 we briefly discuss the implication of the cur- 
rent Fermi results on the detection prospect of PeV neutrinos from 
GRBs. We summarize our results in section 5. 



2 INTERPRETING THE LACK OF GEV EXCESS IN 
MOST GRBS AND THE DELAY OF THE ARRIVAL OF 
THE > 100 MEV PHOTONS 

In the leading internal shock model for th e prompt emission 
dNaravan et al. 1 1 19921; IPaczvnski & Xu 1 1 19941 : iRees & Meszaros 
1 1994 iDaigne & Mochkovitchl 1 19981) . the ultra-relativistic out- 
flows are highly variable. The faster shells ejected at late times 
catch up with the slower ones ejected earlier and then power 
energetic forward/reverse shocks at a radius 7?i nt ~ 5 x 
10 13 (ri/300) 2 (ft v /10 ms) cm, where Yi is the initial Lorentz 
factor of the outflow and St v is the intrinsic variability timescale. 
Part of the shock energy has been used to accelerate elec- 
trons and part has be en given to the magnetic field. If the 
outflow is magnetized dUsov 1 1 19921: iDuncan & Thompson 1 19921 : 
iLvutikov & Blandford l2003l : lGiannios & Spruit l2005h . we call the 
shocks generated in t he collisions within the outflow the m agne- 
tized internal shocks dSpruit et al. 1 1200 ll : iFan et al! 2004b]). The 
synchrotron radiation of the shock-accelerated electrons may peak 
in soft gamma-ray band and then account for the observed prompt 
emission. This model has been widely accepted for the follow- 
ing good reasons: (1) For an ultra-relativistic outflow moving 
with an initial Lorentz factor Li, the velocity (in units of c) is 
ft = y/i - 1/rf . A small velocity dispersion Sft ~ A/(2rf ) 
will yield a very different Lorentz factor. As a result, internal 
shocks within the GRB outflow seem inevitable. (2) In the nu- 
merical simulation of the collapsar launching relativistic outflow, 
people found highly variable energy de position in the polar re- 
gions in a timescale as short as ~ 50 ms (MacF adven & Wooslevl 
Il999l) . (3) This model can naturally account for t he variability that 
is we ll detected in prompt gamma-ray emission dKobavashi et al.l 
1 1997b . On the other hand, this model usually predicts a fast cool- 
ing spectrum Fv oc v^ 1 ^ in the X-ray band. Howe ver, the data 
analy s is finds a typical X-ray spectrum F v oc u° dPreece et al. I 
l200d iBand etaiTl 19931) . Such a divergence between the model 
and the observatio n data is the so-called "th e low energy spec- 
tral index crisis' ' dGhisellini & Celotti I [l999T) . Another potential 
disadvantage of the internal shock model is its low efficiency of 
converting the kinetic en ergy of the outflow into prompt emis- 
sion (e.g.. lKumar]|l999ri . Among the various solutions put for- 
ward, jij3hiuj^ibl£j^c£njffl^^ 

(e.g. IRees & Meszaros Il2005l : IPe'er et al. Il200d : iThompson et al. I 
2007). The idea is that the thermal emission leaking from the 



photosph ere is the dominant component of the prompt sub-MeV 
photons dThompson|[T994l : iMesziros & Reesll200d : lRvdell2005l : 
lloka et aU T2007). The nonthermal high energy emission is likely 
the external inverse Compton (EIC) radiation of the internal shock- 
accelerated e lectrons cooled by the t hermal photons fr o m the 
photosphere dRees & Meszaros 1 120051 : iPe'er et al. 1 120051 120061 : 
IThompson et al.l 2007t) . If the electrons are accelerated by grad- 
ual magnetic energy dissipation rather than by internal shocks, 
it is called the p hotosphere-gradual magnetic dissipation model 
dGiannios|[2007l) . There is an increasing interest in these two 
kinds of models since: (1) In the spectrum analysis people did 
find evidence s for a thermal emission component in dozens of 
bright GRBs dRvde II 20051 : iRvde et al. II200I : iRvde & Pe'erll2009l : 
iMcGlvnn et al. Il2009h . (2) The emission from the photosphere can 
naturally account for the temporal b ehaviors of the temperature 
and flux of these thermal radiation dPe' erll2008l) . (3) The over- 
all sp ectrum of the prompt emission ca n be reasonably interpreted 
(e.g.. IPe'er etai~ll2006l : iGiannios 1120071) . (4) The GRB efficiency 
can be much highe r than that of the internal shock model (see 
IRvde & Pe'er 2009, and the references therein). 

In this section we test these four models with the current Fermi 
GRB data. It is somewhat surprisingly to see that none of these 
models have been ruled out. 



2.1 Explaining the lack of GeV spectrum excess in most 
GRBs 

2.1.1 The standard internal shock model 

In this model, the outflow is baryonic and the thermal emission dur- 
ing the initial acceleration of the outflow is ignorable. The prompt 
emission is powered by energetic internal shocks. There are three 
basic assumptions, (i) e c , eb, e P fractions of shock energy have 
been given to electrons, magnetic field and protons, respectively 
(note that e c + £b + e P = 1). (ii) The energy distribution of the 
shock-accelerated electrons is a single power-law. (iii) The prompt 
soft gamma-ray emission is attributed to the synchrotron radiation 
of the shocked electrons. 

For internal shocks generating at i?i nt , the typical random 
L orentz factor of the e lectrons can be estimated as (see section 4. 1. 1 
of IFan & Piranl d2008h for details) 

~ 760 (l+y ssc ) 1/4 L- 1 n ^ 2 </ 2 13 (l+^) 1/2 ( ep /300keV) 1/2 , 

(1) 

where e p = hv m is the observed peak energy of the synchrotron- 



radiation spectrum (vF v ), h is the Planck's constant, L syn is the 
synchrotron-radiation luminosity of the internal shock emission, 
Y ssc ~ [-1 + A /TT4ee/( | l + fl 2 )e B "l/ 2 is t he regular SSC pa - 
rameter dSari & Esinll200ll : lFan & Pirarj|2008l : IPiran et al. II2009I) . 
and g ~ 7^ m £ p /rjm c c 2 . The SSC in the extreme Klein-Nishina 
regime (g S> 1) is very inefficient. If that happens the non-detection 
of GeV spectrum excess in most Fermi GRBs can be naturally 
explained. With the typical parameters adopted in eq.Q} we have 
g ~ 1, for which the SSC may still be important (i.e., Y ssc ^ 1). In 
this work the convenience Q x = Q/10 x has been adopted in cgs 
units except for some specific notations. 
The SSC radiation will peak at 



hv„ 



1 ' 

27 , 



I + 2.9 



220GeV(l + 2 5 r 1 (l + Y ssc ) 



1/2 



Cn,52-Ri*t,i 3 (l + *)( £p /300 keV) 2 



(2) 



Taking into account the energy loss of the electrons via in- 
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2.1.2 The photosphere— internal shock model 

iThompsonl dl994h proposed the first photosphere model for the 
prompt gamma-ray emission, in which the nonthermal X-ray and 
gamma-ray emission are attributed to the Compton upscattering of 
the thermal emission by the mildly relativistic Alfven turbulence. 
The spectra of GRBs c an be nicely reproduced (see also lPe'er et al. I 
2006: Giannios 120071) . It is, however, difficult to e xplain the energy 
dependence of the width of the gamma-r ay pulse l Fenimore et al. I 
Il995t iNorris et al. Ifl996l) . As shown in iThompson et al. I i20ffm , 
such a puzzle may be solved in the photosphere— internal shocks 
model, in which the sub-MeV emission is dominated by the ther- 
mal emission of the fireball and the nonthermal tail is the EIC ra- 
diation of the electro ns accelerated in internal shocks at a radius 
fiint ~ 10 cm (e.g., Meszaros & Rees 1l2000l : lRees & Meszaros I 
120051 IPe'er et aT71 l2006h . In this model the shock accelerated- 
electrons take a power-law energy distribution dn/dr/' e oc j' e p for 



Figure 1. A possible interpretation of the non-detection of the SSC compo- 
nent by Fermi satellite in the internal shock model. The high energy photons 
with an energy > hu cut have been absorbed by the soft gamma-rays and 
energetic pairs are formed. The pairs will lost their energy through syn- 
chrotron radiation and/or inverse Compton scattering and then produce soft 
gamma-rays. 



verse Compton scattering on prompt soft gamma-rays, the cooling 
Lorentz factor can be roughly estimated as 

7c,c ~ 0.03I/~n :52 i?mt, 13^.5- (3) 

In reality 7 C , C is always larger than 1. The derived 7 C , C < 1 just 
means that the electrons have lost almost all their energies and are 
sub-relativistic. 

Prompt high energy photons above the cut-off frequency v cn t 
will produce pairs by interacting wit h softer photons and wi ll not 
escape from the fi reball. Following iLithwick & Sari and 
iFan & PiraiJ ( l200ot) . we have 



hv c 



2 GeV (1 + *) _1 (e P /300 keV) (2 - p)/p Z,- 2 n /? 2 



5t 



2/p r (2p+8)/p 
-2 1 i,2.5 



(4) 



The SSC radiation spectra can be approximated by F Vsac oc 
i/~ 1/2 for v m < v < f m>S sc, and F„ ssc oc v~ p/2 (,F Vaac oc v~ p ) 
for v > f m ,ssc and g ^ 1 (g ^> 1). The energy ratio of the SSC 
radiation emitted below u cut to the synchrotron radiation in the en- 
ergy range v m < v < max{vM, v cut } can be estimated as 



R 



(p - 2) (v cut I v 

m,ssc) ^ssc 1 



(5) 



where vyi « 30 r\(l + z) MeV is the maximal synchrotro n 
radiation frequency of the shocked electrons JCheng & Wei 1996). 

For Y sac ~ 1, p ~ 2.5 (corresponding to the typical 7— ray 
spectrum F v oc i/ -1,25 for hv > e p ) and v cut <C f m , ssc , we 
have R < 1. Therefore there is no GeV excess in the spectrum, 
in agreement with the data. In other words, the non-detection of a 
significant high energy component is due to a too large hv m ^ sc ~ 
TeV and a relative low hv cut ~ GeV (see FigQ]for a schematic 
plot). The other possibility is that (1 + g 2 )EB > 4e c , for which 
Y S3C ~ 0(1), i.e., the SSC radiation is unimportant and can be 
ignored. 



trons cool down rapidly. As shown in lSari et aT l( fl99i) . in the pres- 
ence of steady injection of electrons, the energy distribution can be 
approximated by iV 7o oc 7'~' p+1 ^ for 7 C > max{7 cc , 7^ m } and 



N- 



7c OC 7 , 



for 7e, c < 7c < 



Under what conditions can 



shock acceleration generate a particle distribution with j' e m ~ 1, 
with a significant fraction of the outflow energy deposited in the 
nonthermal particles? There could be two ways. One is to assume 
that electron/positron pair creation in the outflow is so significant 
that the resulting pairs are much more than the electrons associated 
with the protons. As a result, the fraction of shock energy given to 
each electron/positron will be much smaller than that in the case of 
a pair -free outflow and a 7c, m ~ 1 is achievable (Thompso n et al. I 
2007). The other way is to assume that the particle heating is con- 
tinuous. In the internal shock scenario, this could happen if an out- 
flow shell consists of many sub-shells and the weak interaction be- 
tween these sub-shells may be able to produce multiple shocks that 
can accelerate electrons continually with a very small y' em . 

Since both 7 c m and 7 C C are ~ 1, the EIC spectrum should be 



F„ 



-p/2 



in the MeV-TeV energy range and there is no GeV 



excess forp ~ 2.5, consistent with the Fermi data. 



2. 1.3 The magnetized internal shock model 

If the unsteady GRB outflow carries a moderate/small fraction of 
magnetic field, the collision between the fast and slow parts will 
generate strong internal shocks and then produce energetic soft 
gamma-ray emission. As usual, the ratio between the magnetic en- 
ergy density and the particle energy density is denoted as a. In 
the ideal MHD limit, for a > 1 just a very small fraction of the 
upstream energy can be converted into the downstream thermal en- 
ergy. Therefore the GRB efficiency is very That's why peo- 



1 For the Poynting-flux dominated outflow (i.e., a 3> 1), lUsov 1 jl994l) 

proposed that at the radius r MHD ^6x 10 L 52 cr^ tv. m q-t : o K Cm, 

the MHD condition breaks down and large scale electromagnetic waves 
are generated, where L is the total luminosity of the outflow and t„, m 
is the minimum variabili ty timescale of the central engine (see also 
Zhang & Me szaros 1 12002|) . The particles are accelerated and the syn- 
chrotron radiation of the ultra-relativistic electrons pea ks at f m ~ 5 X 
10 19 CT 3(e c /0.2) 2 [3(p-2)/(p- l)] 2 r 2 .5t l ,,m_ 3 Hz iFan et al. |20o3) . 
provided that a significant part of the magnetic energy has been converted 
to the thermal energy of the particles. In this scenario, the SSC radiation is 
expected to be very weak because in the rest fra me of the e lectrons having 
7 c , m ~ 10 4 (e c /0.2)[3(p - 2)/(p - l)]cr 2 fran et al. iPiobl) . the soft 
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pie concentrate on the internal shock s with a magnetization a ^ 1 
dSpruit et al. ll200ll ; lFan et"ail2004bT) . 

For the magnetized internal shocks with a 0.1, no signif- 
icant high energy emission is expected since: (a) The SSC emis- 
sion of the internal shocks is weak. Therefore there is no distinct 
GeV excess in the spectrum, (b) The synchrotron-radiation spec- 
trum may be very soft, which renders the detection of GeV pho- 
tons from GRBs more difficult. The reason is the following. For 
an isotropic diffusion and a relativistic shock, the electron energ y 
distribution index can be estimated by (K eshet & Waxman 1120051) 



(3/3 u 



2/Juifld + 0d)/O3a - #0 



(6) 



However, in the presence of a large scale coherent magnetic 
field, the diffusion i s highly anisotropic rather than isotropic 
(Morli ni et al. I l2007ah . There are thus corrections to eq.(|6). But 
as long as the scattering is not very forwa rd- or backward-peaked, 
these corrections are small ( Keshet jgoodi Taking into account the 
anisotropic correction, Morlin o et al. T l2007bl) found a spectrum 
steep to p ~ 3 for a ~ 0.05. In the ion-electron shock simulation, 
the acceleration of particles at the un-magne tized shock front is a 
lot m ore efficient than that with a a ~ 0.1 ( Sironi & Spitkovsky 
2009). Motivated by these two possible evidences, we adopt eq.{6} 
to estimate the spectral slope of accelerated particles at the magne- 
tized shock fronts. The validity of our approach can be tested by 
the advanced numerical simulations in the future. 

In the case of a — 0, for an ultra-relativistic shock, /3u 1 
and Pd — * 1/3, we have p — > 2.2 2. But for an ultra -relativistic 
magnetized shock, j3 u — > 1 and (e.g jFan eta l. 2004b) 



3 d 



^(i + x + \A + i4x 



■ X s 



(7) 



where x = + a )> please note that a is measured in the up- 
stream. Note in this work we just discuss the ideal MHD limit, 
i.e., there is no magnetic energy dissipation at the shock front. For 
< a < 1, we have p ~ (4.22 - 2a) /{I - 2a) - 2 > 2.22. For 
a > 1, we have f3 d -> 1 - 1/2ct and p ~ 4a - 1 » 2.22. Corre- 
spondingly, the electron spectrum is very soft or even thermal-like. 
Adopting /3 U ~ 1 and substituting a ~ (1, 0.5, 0.1, 0.01) into 
eqs.l |6|7K we have p ~ (6.6, 4.5, 2.7, 2.3). The (ve ry) soft high 
energy spectra of some GRBs (e.g. JPreece et al. l200d see also our 
Fig|2]for the Fermi GRBs) may be interpreted in this way. 

If the weak prompt high energy emission of GRBs is indeed 
attributed to the magnetization of the outflow, one can expect that 
the smaller the p, the stronger the high energy emission. The ongo- 
ing analysis of the LAT data will test such a correlation. 



2.1.4 The photosphere— gradual magnetic dissipation model 

iGianniosI J2007I) calculated the emission of a Poynting-flux- 
dominated GRB outflow with gradual magnetic energy dissipation 
(reconnection). In his scenario, the energy of the radiating electrons 
is determined by heating and cooling balance. The mildly relativis- 
tic electrons stay thermal throughout the dissipation region because 
of Coulomb collisions (Thomson thick part of the flow) and ex- 
chang e of synchrotron photons (Thomson thin part). Rather simi- 
lar to lThompson I d 1994b . the resulting spectrum naturally explains 



gamma-rays have an energy 7^ m /ii^ m /r; 2> m e c 2 , i.e., the SSC is in the 
extreme Klein-Nishina regime and is very inefficient. The non-detection of 
high energy emission from most GRBs is, of course, consistent with this 
model. 




-1.5 -2.0 -2.5 -3.0 -3.5 

high energy power-law index 

Figure 2. The distribution of the high energy power-law index (/?Band) of 
GRBs detected by Fermi satellite from 10 August 2008 to 31 March 2009. 
The data are taken from http://www.batse.msfc.nasa.gov/gbm/circulars/ 
(see also http://gcn.gsfc.nasa.gov/gcn3— archive. html I and are preliminary. 
The GRBs having a single power-law spectrum in the GBM energy range 
are excluded in the statistics. Following Pre ece et al. I j200Ql) we take 
/^Band = — 4 for all bursts that can be better fitted by a power-law 
function with an exponent ial high energy cutoff rather than by the Band 
function feand et alj|l993l) . We find that a goo d fraction of bursts have a 
p ~ -2( / 3 Ban d + 1) > 2.22, consistent with lPreece et al. I l2000h . These 
soft spectra may be attributed to the magnetized internal shocks. 



the observed sub-MeV break of the GRB emission and the spec- 
tral slopes. In this scenario, different from the magnetized internal 
shock mod el, the higher the initial a, the harder the spectrum (see 
the Fig.2 of|G tanntos 1 J2007b for illustration). For an initial a ^ 40 
(corresponding to the baryon loading L/Mc 2 ~ a 3//2 ^ 250, 
where M is the mass loadi ng rate), the resulting > 10 MeV spec- 
trum is very soft (see also Drenkhahn 1 l2002l : lDrenkhahn & Spruitl 
2002), accounting for the failed detection of the GeV spectrum ex- 
cess in most GRBs. 



2.2 Interpreting the delay of the arrival of the > 100 MeV 
photons 

In both the collap s ar an d the compact sta r merger models for 
GRB s (see IPiranl 1 19991 : iMeszaros I 120021 : Izhang & Meszaros I 
120041 for reviews), the early outflow may suffer more serious 
baryon pollution and thus have a sma ller Ti than the late ejecta 
( Zhang, Wooslev & MacFadven 1 120041) . This may explain the de- 
lay of the arrival of the > 100 MeV emission since as long as 



A sC Ti.c = 180 (l+*) 2 *+ B ( 



hv cu 



100 MeV 



) 2p + S ( 



300keV y 



p-2 1 



syn,51 



the > 100 MeV photons can not escape from the emitting region 
freely and thus can not be detected. 

In the photosphere-gradual magnetic dissipation model, 
a small Ti implies a low initial magnetization of the out- 
flow, for which the high energy spectrum can be very soft 
jDrenkhahn & Spruit 1 120021 : IGiannios I 120071) . In the magnetized 
internal shock model, the delay of the onset of the LAT observa- 
tion indicates a larger magnetization of the early internal shocks if 
Ti > r . 

In the collapsar scenario, before the breakout, the initial out- 
flow is choked by the envelope material of the massive star (Zhang 
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et al. 2004). The ultra-relativistic reverse shock may be able to 
smooth out the velocity/energy-density dispersion of the initial 
ejecta. So the internal shocks generated within the early/breakout 
outflow may be too weak to produce a significant non-thermal ra- 
diation component. The early emission is then dominated by the 
thermal component from the photosphere and may last a few sec- 
onds (provided that the chocked material has a width comparable 
to that of the envelope of the progenitor). The outflow launched af- 
ter the breakout of the early ejecta can escape from the progenitor 
freely and the consequent internal shocks can be strong enough to 
produce energetic non-thermal radiation. The photosphere-internal 
shock model therefore might be able to naturally account for the 
delay in the onset of the LAT observation. 

In summary, before and after the onset of the > 100 MeV 
emission, it seems the physical properties of the outflow have 
changed. 



z' (W 




3 THE LINEAR POLARIZATION SIGNAL OF THE 
PROMPT 7 RAY EMISSION 

As discussed in section 2, the failed detection of the GeV spectrum 
excess in most GRBs can be understood in either the standard in- 
ternal shock model or several alternatives. Therefore we need inde- 
pendent probes to distinguish between these scenarios. Our current 
purpose is to see whether the polarimetry in gamma-ray band can 
achieve such a goal. In this section, we firstly investigate the lin- 
ear polarization property of the photosphere-internal shock model 
(the results may apply to the photosphere-gradual magnetic dissi- 
pation model as well) since it has not been reported by others yeQ 
We then briefly discuss the linear polarization signals expected in 
the magnetized internal shock model and in the standard internal 
shock model since they have been extensively discussed in the lit- 
erature (e.g iLvutikov et al. 1120031: i Granot 1120031; IWaxman 1120031 ; 
iNakar et al. ll2003l : lFan et al. booalToma et al. Il2009h . 



3.1 Linear polarization signal of the photosphere-internal 
shock model 

In this work we assume an uniform outflow. At any point in the 
outflow there is a preferred direction, the radial direction, in which 
the fluid moves. We choose the Z'-direction of the fluid local frame 
coordinate to be in that direction. The Y'— direction is chosen to 
be within the place containing the line of sight (i.e., the scattered 
photon k' ) and the Z'— axis (see Fig[3]>. In this frame, the incident 
photons (k') are along the Z'— direction. 

As usual, we assume that the electrons are isotropic in the co- 
moving frame of the emitting regiorff]. The incident photons are 
the thermal emission from the photosphere and are unpolarized. 



2 In the final stage of preparing the manuscript, the author was informed 
by K. Toma and X. F. Wu that the POET (Polarimeters for Energetic 
Transients) group discussed the polarization property of the photosphere- 
internal shock model independently. The preliminary results were included 
in their white paper. 

3 In the current scenario, the electrons are heated by the internal shocks 
and an isotropic distribution in the rest frame of the emitting region may 
be a good approximation. Such a distribution, however, will be modified by 
the cooling of the electrons via EIC scattering on the thermal photons from 
the inner part of the outflow. Detailed numerical simulation is needed to see 
such a modification and then its influence on the inverse Compton radiation, 
which is beyond the scope of this work. 



Figure 3. The coordinates in the comoving frame of the emitting region. 
The polarization vector El' is along the X'— direction for a positive P given 
in eq.(TT} otherwise it would be along X' X k' Q . 



The energy of the incident and t he scattered photons a r e hv' s e and 
hv' , respectively. As sh own in [Aharonian & Atovaril dl98ll) and 
iBerestetskii et al. I Jl982h : (I) In the cloud of the isotropic electrons 
of energies ■y' e m e c 2 , the spectrum of photons, upscattered at the 
angle 9' relative to the direction of the seed photon beam, can be 
approximated by 



dN-, 



x {AAl-AA + B ),{&) 



dtdv'dQ.' 167T7' 2 v' se 2%Q 

where $ = (1 - l/ilT 1/2 and 

a - 1 / 1 1 \ 

2dy e x Ti 7 2 



So 



% 71 
71 + T 2 ' 



_ , 1 + x + x5cos9' — Sx 2 
71 = (1 — cosO) 



T 2 = (1 - cos 



Q 2 

1 + x — x8 cos 6' + 5 
' Q 2 ' 



Q = \/l + x 2 - 2x cos 6", 

8 = hv' sc /(~/ c m c c 2 ), x = v'/v' se . (9) 
Please note that x ranges from Xm to x M that are given by 



= 1 



{A T I'J'cVYci 1 + S) 2 ^ ~ cos^') 2 + sin 2 6'} 



l + 2<5 7 ' e 2 (l 



') +l' 2 J 2 {l - cos6") 2 



where A= (1 - cos 0')Ye [0 e - 5 - <5 2 (1 - cosd')]. 
(II) The polarization degree is 

B - 4A + AA 2 ■ y ' 

For the relativistic electrons (i.e., /3 e — > 1), eq.{8]( and eq.dl lb 
take the simplified forms 

dN-y 



dtdv'dQ,' 



3a T c n v i du' sc £ 2 
[! + 



2£ 



16^7'. 
2C 2 



b 2 (l-0 2 > 



2(1-0 be(l-0 

(12) 
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Figure 4. Sketch of the geometrical set-up used to compute the polarization 
signal. We take the L.o.S as the z— axis. The y-axis (a;— axis) is within 
(perpendicular to) the plane containing the line of sight and central axis of 
the ejecta. 



2« 



MW) 



2C 



1 + 



2? 



2(i-« Mi-e) 



2£ 2 
+ ^(i-€) 2 



(13) 



where £ = his' /(j c m c c 2 ), bg = 2(1 — cos 7' e hv' sl , / {m^c 2 ), 
and /u4 < W < 7' c m c c 2 6 9 /(l + b e ). 

Since the emitting region is moving relativistically, the angle 
8' (in Figf3]l corresponding to the line of sight (L.o.S) is given by 

ros0' = (cos0-A)/(l-A«Mi0), (14) 

where is the angle between the line of sight and the emitting point 
(measured in the observer's frame). 

The azimuthal angle <j> varying from to 2n is defined in 
FigE] The polar angle 9 ranges from to 9 V + 8j . The angle be- 
tween the vector (sin 6 cos <j>, sin 9 sin cj>, cos 9) and the central 
axis of the ejecta (C.A. in FigQ is denoted as G and is given by 

cos — ~ sin 8 V sin 9 sin (j> + cos 9 V cos 9. (15) 

Please bear in mind that in the following radiation calculation, the 
flux is set to be zero if cos O < cos 9j because these points (9, <j>) 
are outside of the cone of the ejecta. 

The EIC radiation flux in the observer frame is 



■Fi/.Eic oc / V 3 hv' 



dtdv'dn-^'^ (16) 

where v = 2V/(1 + z), V = \Ti{\ - $ cosS)]" 1 is the Doppler 
factor, and Q. is the solid angle satisfying dQ, — si 
The polarized radiation flux is 



Qvf eic oc / D z hv'P cos 20- 



C2./V.V 



-N^d-y'dQ. (17) 



The polarization degree of the EIC emission is 

Pt/,mc = IQ^.Eicl/i^.Eic- (18) 

One can see that a non-zero net polarization is expected as 
long as 9 V > 0. In the numerical example, we assume that the seed 
photons have a thermal spectrum (as suggested in the photosphere 
model) 



hv^ B /kT> 



(19) 



10" 10'" 10" 10" 

frequency (Hz) 

Figure 5. The degree of the linear polarization of the EIC component as 
a function of frequency. The half-opening angle of the ejecta 8j and the 
bulk Lorentz factor T; are assumed to be 0.1 and 300, respectively. The 
viewing angle 8 V is marked in the plot. The solid and dashed lines are for 



the electron distribution N. 



70 oc 7o 



' and N, 



oc 7^ 



, respectively. 




10'" 10" 
frequency (Hz) 



Figure 6. The spectrum of one simulated GRB (the upper panel) and the 
expected linear polarization degree of the total emission as a function of the 
photon energy (the lower panel). In the calculation we take 8 V = 0.102. 
Other parameters are the same as those used in Fig|6]in the case of jV-y e oc 



where kT' fa kT/2Ti is the temperature (measured in the rest 
frame of the emitting region) of the thermal emission. In the cal- 
culation we take Fi ~ 300 and kT ~ 100 keV. The electron 



distribution is taken as 7V 7c oc 7' 



-a+p) 



oc 7 , 



for 7 ' c > 2, 



otherwise 7V 7o = 0. For comparison purpose we also consider the 
case of iV 7o oc 7'~ 2 . The numerical results are presented in Fig(5] 
One can see that the polarization degrees expected in these two 
representative cases are only slightly different. We also find that a 
moderate linear polarization level (P„,eic > 10%) is achievable 
only for 9 V > 9 S + 1/(3^). 

Currently the prompt emission consists of a thermal and a non- 
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thermal components. The thermal component with the flux F v ^ is 
expected to be unpolarized while the nonthermal EIC component 
may have a high linear polarization level. The observed polarization 
degree 



Table 1. Distinguish between the models with the polarimetry data. 



\Q».\ 



F V|t h 



should be strongly frequency-dependent. Roughly speaking, the 
linear polarization degree is anti-correlated with the weight of the 
thermal component. With an energy ~ kT, the emission is domi- 
nated by the thermal component and PkT.obs is low. For hl> 2> kT, 
the emission is dominated by the EIC component and P„,ot>s ~ 
Pv.eic, as illustrated in Fig(6] This unique behaviojf] can help us to 
distinguish it from other models. 

The probability of detecting a moderate/high linear polariza- 
tion degree (TZ po i), however, is not high (Please note that we do 
not take into account the weak events for which a reliable po- 
larimetry is impossible). On the one hand, a high linear polarization 
level is achievable only for 8 V ^ 8j + 1/3IV On the other hand, 
8 V — 9j < 1/Tj is needed otherwise the burst will be too weak to 
perform the gamma-ray polarimetry. For Fi#j 3> 1, we have 

ftpoi ~ 4/(31^) « 5% rr^tfr^. (20) 

During the revision of this work, McGlyn n et al. I d2009h re- 
ported their analysis on the spectrum and the polarization proper- 
ties of GRB 061122. They found out that the spectrum was better 
fitted by the superposition of a thermal and a non-thermal com- 
ponents and the photons in the "thermal" emission dominated en- 
ergy range had a (much) lower polarization level than those in the 
higher energy band. These two characters are in agreement with 
the photosphere-internal shock model (or the photosphere-gradual 
magnetic dissipation model). 



3.2 Linear polarization level expected in the standard 
internal shock model 

In the standard internal shock model, the polarization of the syn- 
chrotron radiation depends on both the poorly known configuration 
of magnetic field generated in internal shocks and the geometry 
of the visible emitting region. Assuming a r andom magnetic field 
that remains planar in the plane of the shock. IWaxman U2003I) and 
Nakar^etaT] d2003h showed that a high linear polarization level can 
be obtained when a narrow jet is observed from the edge, like in the 
photosphere-internal shock model. For the jets on-axis (9 V ^ 8j), 
the linear polarization degree is low (see also Gruzinov 1999 and 
Toma et al. 2009). The detection probability of a moderate/high 
linear polarization degree can also be estimated by eq.J20b. 



3.3 High linear polarization degree expected in the 
magnetized internal shock model 

For the magnetized internal shock model, the prompt soft 7— ray 
emission is attributed to the synchrotron radiation of the electrons 
in ordered magnetic field and a high linear polarization level 



4 For the synchrotron radiation of the electrons moving in a random mag- 
netic field (the standard internal shock model) or in an ordered magnetic 
filed (e.g., the magnetized internal shock model), before and after the peak 
of the spectrum, the polarization degree changes because the polarization 
properties depend on the p rofile of the sp ectrum. However, such a depen- 
dence is weak, as shown in Gr anot 1 120031) . 



model 


unique polarization property 




standard internal shocks 




< 10% 


photosphere-internal shocks^ 


strongly frequency-dependent^ 


< 10% 


magnetized internal shocks 




~ 100% 



tin the photosphere-gradual magnetic dissipation model, very similar 
polarization properties are expected. 

''As shown in Fig|6] the polarization degree is anti-correlated with the 
weight of the photosphere/thermal component. 



is expected jLvutikovetaLll2003l ; iGranot I l2003h . The physical 
reason is the following. The magnetic fields from the central 
engine are likely frozen in the expanding shells. The toroidal 
magnetic field component decreases as while the poloidal 

magnetic field component decreases as R~ 2 . At the radius of the 
"internal" energy dissipation (or the reverse shock emission), the 
frozen-in field is dominated by the toroidal component. For an 
ultra-relativistic outflow, due to the relativistic beaming effect, 
only the radiation from a very narrow cone (with the half-opening 
angle ^ 1/Ti) around the line of sight can be detected. As long as 
the line of sight is off the symmetric axis of the toroidal magnetic 
field, the orientation of the viewed magnetic field is nearly the 
same within the field of view. The synchrotron emission from such 
an ordered magnetic field therefore has a preferred polarization 
orientation (i.e. perpendicular to the direction of the toroidal field 
and the line of sight). Consequently, the linear polarization of the 
synchrotron emission of each electrons could not be effectively 
averaged out and the net emission should be highly polarized. 
The detection prospect of a high linear polarization degree is very 
promising (i.e., TZ po \ ~ 100%). The above argument applies to 
the reverse shock emission as well if the outflow is magnetized 
dFan et al. Il2004ij) . 

As summarized in TabQ] for the magnetized internal shock 
model, a high linear polarization level should be typical, while for 
two other models a moderate/high linear polarization degree is still 
possible but much less frequent. So the statistical analysis of the 
GRB polarimetry results may be able to distingui sh the magnetize d 
internal shock model from the others (see aslo lToma et al. 1 2009). 
In the photosphere-internal shock model the polarization degree is 
expected to be strongly frequency-dependent. Such a remarkable 
behavior, if detected, labels its physical origin. 

Indeed there were some claims of the detection of high 
linear p olarization degree in the soft 7— ray emission of GRB 
021206 dCoburn & Boggsl 120031 how ever see Rutledge & Fox 
2004), G RB 930131, GRB 960 9 24 dwillis et al. I |2005t) . GRB 
041219 A jMcGlvnn et all 120071 : iGotz et al. I I2OO9J) , and GRB 
061122 dMcGlvnnetal.1120091) . These results are consistent with 
each other as the errors are very large. The situation is inconclusive 
and additional data is needed to test these results. Measuring polar- 
ization is of growing interest in high energy astronomy. New tech- 
nologies are being invented, and several polarimeter projects are 
proposed, such as, in the ga mma-ray band, ther e are the Advanced 
Comp ton Telescope Mission dBoggs et al.l2006f) . POET faill et al. I 
2008) and others (see Tom a et al. 1120091 for a summary). So in the 
next decade reliable polarimetry of GRBs in gamma-ray band may 
be realized and we can impose tight constraint on the models. At 
present, the most reliable polarimetry is in UV/optical band (e.g., 
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ICovino et al. Il999l ; lwiiersetal.ll999h . The optical polarimetry of 
the prompt emission and the reverse shock emission require a quick 
response of the telescope to the GRB alert. This is very challeng- 
ing. Mundell et al. (2007) reported the optical polarization of the 
afterglow, at 203 sec after the initial burst of 7— rays from GRB 
060418, using a ring polarimeter on the robotic Liverpool Tele- 
scope. Their robust (90% confidence level) upper limit on the per- 
centage of polarization, less than 8%, coincides with the fireball 
deceleration time at the onset of the afterglow. Such a null detec- 
tion is, however, not a surprise because for this particular burst the 
reverse shock emissio n is too weak to o utshine the unpolarized for- 
ward shock emission dJin &Fanll2007l) . Quite recently, the robotic 
Liverpool Telescope performed the polarimetry measurement of the 
reverse shock emission of GRB 090102 (Ko b ayashi 2009, private 
communication). Following IFan et al. I j2002h . lZhang et"aI1 d2003l) 
and lKumar & Panaitescu 1 1 20031) . it is straightforward to show that 
the reverse shock of GRB 090102 is magnetized. Consequently the 
optical flash is expected to be highly polarized. If confirmed in 
the ongoing data analysis, the magnetized outflow model for some 
GRBs will be favored. 



4 IMPLICATION ON THE DETECTION PROSPECT OF 
PEV NEUTRINO EMISSION 

The site of the prompt 7— ray emission may be an ideal place 
accel erating protons to ultra-high energy dVietrilll995l : I Waxman 
[l99l) . These energetic protons can produce high-energy neu- 
trinos via photomeson inte raction, mainly through A-resonance 
(Wax man & Bahcail"| [T997). The resulting neutrinos have a typi- 
cal energy £„, oba ~5x 10 14 eV r|. 6 [(l + z) 2 e 7 ,ob s /l MeV]" 1 . 
Significant detections are expected if GR Bs are the main s ource 
of ultra-high energy cosmic rays dWaxman" &Bahcall|[l997l) . The 
underlying assumption is that the proton spectrum is not signifi- 
cantly softer than dN/dE oc E~ 2 . The current Fermi observa- 
tions do not provide an observational evidence for such a flat par- 
ticle spectrum. Below we discuss the detection prospect of PeV 
neutrinos implicated by the non-detection of high energy emis- 
sion from most GRBs. In the magnetized outflow model, the ac- 
celeration of a significant part of protons to energies > 10 16 eV 
is highly questionable because of the resulting soft proton spec- 
trum. In the photospher e-internal shock model, 7^ m ~ 1 is needed 
dThompson et al. ||2007|) . The efficiency of accelerating protons to 
very high energy depends on the mechanism of the particle heat- 
ing. For example, in the case of multiple internal shocks, each pair 
of internal shocks are expected to be very weak since r s h — 1 ~ 
0.04( 7 ; m /5)(e c /0.2)- 1 [3(p - 2)/{p - l)]" 1 , where T sh is the 
Lorentz factor representing the strength of the shock (r s h ~ 1 for 
Newtonian shocks). So the acceleration of the protons to ultra-high 
energy is less efficient than the standard internal shocks. This is 
particularly the case if the acceleration is mainly via second-order 
Fermi process, in which the acceleration of particles depends on 
the shock velocity sensitively. 

Even in the standard internal shock model, the generation of 
10 20 eV protons and the production of PeV-EeV neutrinos may be 
not as promising as that claimed in most literature adopting a pro- 
ton spectrum dN/dE oc E~ 2 . Such a flat spectrum is predicted for 
Newtonian shocks and has been confirmed by the supernova rem- 
nant observations. B ut the typical MeV spectrum F v oc v~ 1,25 
dPreeceetal.ll2000h of GRBs suggests dN/dE oc S" 2 ' 5 , sup- 
posing the accelerated protons and electrons have the same spec- 
trum. The non-detection of > 100 MeV photon emission from most 



GRBs implies soft electron (possibly) and proton spectra. Given a 
proton spectrum dN/dE oc E~ 2 ' 22 that is predicted in the rel- 
ativistic shock acceleration model (the First-order Fermi mecha- 
nism), the kinetic energy of the ejecta needs to be ~ 100 times of 
the 7 — ray radiation energy if GRBs are in deed the main s ource of 
the observed ~ 10 20 eV cosmic rays (see iDermer 1120081 and the 
reference therein). In other words, the GRB efficiency should be 
as low as ~ 1%. If correct, the number of protons at E ~ 10 1 6 
eV will be quite a few times what assumed in lGuetta et al. I d2004l) . 
Correspondingly the PeV neutrino flux will be higher. However, 
the current afterglow modeling usually yields a typical GRB ef- 
ficiency - 10% (e.g., iFan & Piranl 120061; IZhang et al. 1120071) or 
larger (e.g., IPanaitescu & Kumar II2002I : iGranot et al. 112006^ Be- 
low we discuss a new possibility-The proton spectrum is curved. 
In the "low-energy" part, the spectrum may be steepened signif- 
icantly by the leakage of the ver y high energy cosmic rays from 
the ejecta (see lHillas et al. 1 12005. and the references therein). The 
"high-energy" spectrum part may be a lot flatter. For example, in 
the numerical simulation of cosmic rays accelerated in some super- 
nova remnants, a spec trum dN/dE oc E~ 17 at the high ene rgy 
part is obtained (e.g., IVolketal.ll2002l ; iBerezhko et alj|2003h . If 
holding for GRBs as well and GRBs are the main source of the 
10 20 eV cosmic rays, the PeV neutrino s pectrum will be harder 
than that predicted in iGuetta et al. I d2004h . For instance, the neu- 
tron spectra etdN/de v oc el and e 2 dN '/ds v oc e„ 2 in their Fig.3 
will be hardened by a factor o f e° 3 . But the total fl ux may be just 
~ 10% times that predicted in lGuetta et al. I d2004h because in this 
scenario the protons are not as many as that suggested in a flat spec- 
trum dN/dE oc E-' 2 for E < 10 20 eV. 

There is a process, ignored in some previous works, that can 
enhance the detection prospect a little bit. After the pions (muons) 
are generated, the high energy pions (muons) will lose energy via 
synchrotron radiation before decaying, thus reducing the energy of 
the decay neutrinos (e.g. Guetta et al. 2004). As a result, above 



1/2 -l/2 r -l/2 p4 



7,52 



5 <5tv,-2 eV (£% u 



S M /10), 



the slope of the corresponding neutrino spectrum s teepens by 2, 
where L 7 is the luminosity of the 7— ray emission dGuetta et al. I 
2004). However we do not suggest a smooth spectral transition 
around e% or e% ^ because the cooling of pions (muons) will 
cause a pile of particles at these energies (see also Murase & Na- 
gataki 2006). A simple estimate suggests that the number of neu- 
trinos in the energy range (0.5, should be enhanced by a 
factor of ~ 3. A schematic plot of the muon neutrino spectrum in 
the standard internal shock model is shown in Fig|7] 



5 CONCLUSION 

In the pre-Fermi era, it is widely expected that significant GeV 
emission will be detected in a good fraction of brigh t GRBs if they 
are powered by un-magnetized internal shocks (e.g., Pilla & Loeb I 



1 1998 1 ; IPe'er & Waxman 1120041 ; iGupta & Zhang1l2007l ; IFan & Piranl 
2008). The detection of a distinct excess at GeV-TeV energies, the 
SSC radiation component of such shocks, will be a crucial ev- 
idence for the standard fireball model. The no n-detection of the 
GeV spectrum excess in almost all Fermi bursts (Abd oetalj2009T) 
is a surprise but does not impose a tight constraint on the mod- 
els. For example, in the standard internal shock model, the non- 
detection can be attributed to a too large /i^ m ,ssc ~ TeV and a rela- 
tive low hv cu t ~ GeV. Some alternatives, such as the photosphere- 
internal shock model, the magnetized internal shock model and the 
photosphere-gradual magnetic dissipation model, can be in agree- 
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Figure 7. The schematic plot of the PeV muon neutrino spectrum in the 
standard internal shock model. 



Table 2. The physical reasons for the lack of GeV spectrum excess in most 
GRBs. 



model 


the physical reason 


standard internal shocks 


h(u m ,ssc, i/cut) ~ (TeV, GeV) 


photosphere-internal shocks 


very small 7^ m and 7^ c 


magnetized internal shocks 


soft synchrotron spectrum and weak SSC 


photosphere-gradual 




magnetic dissipation 


very small 7^ m and 7^ 



ment with the data, too (see Tab[2]for a summary). We attribute the 
delay in the onset of LAT detection in quite a few Fermi bursts to 
the unfavorable condition for GeV emission of the early outflow 
(see section l2~2l for details). 

With the polarimetry of GRBs people can potentially distin- 
guish between some prompt emission models (see Tab[T]for a sum- 
mary; see also Toma et al. 2009). We show in section [3~Tl that in the 
photosphere-internal shock model the linear polarization degree is 
roughly anti-correlated with the weight of the thermal component 
and will be highly frequency-dependent. Such a unique behavior, 
if detected, labels its physical origin. However, a moderate/high 
linear polarization level is expected only when the line of sight 
is outside of the cone of the ejecta (i.e., 8 V > 9j). In addition, 
f? v — 9} < 1/Ti is needed otherwise the burst will be too weak 
to perform the gamma-ray polarimetry. Consequently the detection 
prospect is not very promising. 

In this work we have also briefly discussed the detection 
prospect of prompt PeV neutrinos from GRBs. The roles of the 
intrinsic spectrum of the protons and the cooling of pions (muons) 
have been outlined. The latter always increases the neutrino num- 
bers at the energies el or e% iVe by a factor of 3. The former, 
however, is uncertain. If the protons have an intrinsic spectrum 
dN/dE oc £?~ 2 - 22 and have a total energy about tens times 
that emitted in gamma-rays, the detection prospect would be a s 
good as, or even better than that presented in I Guetta et al. I J2004I) . 
If the proton spectrum traces that of the electrons, i.e., typically 
dN/dE oc E~ 2 ' 5 , the detection prospect would be discouraging. 
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